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Abstract 

We provide evidence that general relativity is the unique spatially covariant effective 
field theory of the transverse, traceless graviton degrees of freedom. The Lorentz co- 
variance of general relativity, having not been assumed in our analysis, is thus plausibly 
interpreted as an accidental or emergent symmetry of the gravitational sector. 



1 Introduction 



Lorentz covariance is a central pillar of the modern field-theoretic interpretation of general 
relativity (GR). From this point of view, GR is no more and no less than the unique Lorentz 
covariant theory of an interacting massless spin-2 particle [1-3]. In this paper, we show 
that GR can be derived without assuming Lorentz covariance. Our approach relies on the 
weaker assumption of spatial covariance, within the context of the effective field theory of the 
transverse, traceless graviton degrees of freedom developed in [4]. 

In canonical form, GR is a theory of a spatial metric hij subject to first class 1 constraints 
Hp. Each constraint both generates a space-time gauge symmetry and eliminates a physical 
degree of freedom. The Hamiltonian constraint "Ho generates local time reparameterizations 
and eliminates the scalar polarization of the graviton; the momentum constraints % generate 
spatial diffeomorphisms and eliminate the longitudinal polarizations of the graviton. In 3 + 1 
dimensions, the spatial metric has six components subject to four constraints H^, so the 
graviton of general relativity has 6 — 4 = 2 transverse, traceless polarizations. 2 

In this framework, Lorentz covariance arises because the constraints obey the Dirac 
algebra [5, 6], the algebra of the deformations of a spacelike hypersurface embedded in a 
Lorentzian space-time manifold [7]. The project of relaxing the assumption of Lorentz co- 
variance without introducing new gravitational degrees of freedom has been the subject of 
much research, and there are many directions one can take. One can, for example, replace the 
Hamiltonian constraint with a restriction on the form of the Lagrangian [6]. In the ultralocal 
truncation of GR, the Hamiltonian constraint still removes the scalar graviton polarization, 
but the constraints no longer satisfy the Dirac algebra [8] . In the theory of Shape Dynamics, 
the Hamiltonian constraint is replaced by a constraint that generates volume preserving con- 
formal transformations [9, 10]. In the generally covariant version of Hofava-Lifshitz gravity, 
additional fields and constraints are added to the theory to obtain the desired number of 
degrees of freedom [11]. 

Our method for relaxing the Hamiltonian constraint without introducing new degrees 
of freedom is to reduce the number of fields in the theory. Our approach is inspired by a 
technique of conformal geometry, specifically the equivalence class construction underlying 
the notion of a conformal metric. Any two metrics hij and h'^ belong to the same conformal 
equivalence class if there exists a scalar function Q such that h'^ = e n ■ h^. Each conformal 
class defines a geometry up to local changes of scale, and thereby defines a conformal metric. 

^irst class constraints are constraints that all commute with each other. 

2 See [4] for a comprehensive analysis of general relativity as a constrained field theory. 
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Each conformal class also serves to define a unique unit- determinant metric. To see 
this, consider an arbitrary representative metric hij drawn from a conformal equivalence class. 
Following the convention of [8], in d spatial dimensions we define 

hij^e^hij, n = ^\ogh. (1) 

By definition, hij and hij belong to the same conformal equivalence class, and the metric hij has 
unit determinant. Moreover, the metric so defined depends only on the choice of conformal 
class, not on the choice of representative within the class. We have therefore established a 
one-to-one correspondence between conformal metrics and unit-determinant metrics. 

Instead of thinking about gravity in terms of a spatial metric hij subject to a system 
of constraints that kills its conformal mode, we will think of gravity in terms of a unit- 
determinant spatial metric h^ that already lacks an independent conformal mode. This allows 
us to remove the scalar polarization of the graviton while remaining agnostic about precisely 
how and why this polarization is absent. Since it has one fewer component to begin with, a 
spatially covariant unit-determinant metric h^ has the same number of degrees of freedom as 
a space-time covariant metric hij. For example, in 3 + 1 dimensions a unit-determinant spatial 
metric hij has five components, so subjecting such a metric to three momentum constraints 
Hi yields a graviton with 5 — 3 = 2 transverse, traceless polarizations. Insofar as we consider 
metrics with unit conformal factor, our proposal is similar in spirit to unimodular gravity [12] . 

By construction, spatially covariant theories of a unit-determinant metric describe the 
same transverse, traceless graviton polarizations as GR. Though the kinematical state space 
is essentially the same, in principle the dynamical evolution of the graviton degrees of freedom 
could differ dramatically. However, as we will see, demanding a consistent algebra and evolu- 
tion for the momentum constraints singles out GR as the unique spatially covariant effective 
field theory of a transverse, traceless graviton. 

The proof is quite technical in detail, but its outline is simple and can easily be summa- 
rized in a few steps. Our starting point, described in Sec. 2, is the general canonical action in 
d spatial dimensions 

S = J dtd d x (hijjp* -7TH- N*iu) , (2) 

describing a unit-determinant metric hy and its traceless conjugate momentum 7f lJ . The 
scalar function tt^ is the physical Hamiltonian density. Spatial covariance is enforced by the 
Hi momentum constraints, whose general form includes a second scalar function ttk' 

Hi EE -2hijV k f jk - Vi7l K . (3) 
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Here V denotes the covariant derivative with respect to hy. Our set-up therefore depends 
on two functions, the physical Hamiltonian density tth and the momentum constraint density 
7Tk, which a priori are arbitrary functions of t, the phase space variables (7iy,7f y ), and the 
spatial gradient operator <9;. In this language, GR corresponds to a particular form for n H 
and ttk, namely: 

tt£ r = -n(*)7r n , vr£ R = ^ , (4) 



where 



n n = ±^/d{d - I) ^ (f)(2) - l&d-mt) + 2Ae dn ^ . (5) 



Here 0(2) = hikhjtTTijTT , A is the cosmological constant, and Q(t) is a positive-definite, 
monotonic function of time. For our general class of gravity theories, our only simplifying 
assumption about the form of tth and ttk is that they depend on spatial gradients only 
through the Ricci scalar R. 

The allowed form for kh and ttk is constrained by two considerations. First, in order 
to generate a consistent gauge symmetry, the momentum constraints Hi must be first-class 
under the action of the Poisson bracket, 3 

{Hi(x),H a {y)}~0. (7) 

where ~ denotes weak equality, i.e., equality after the imposition of all constraints. The con- 
dition (7) means that the algebra of the constraints is closed, in the sense that the Poisson 
bracket of any two constraints is a linear combination of the constraints. The closure con- 
dition (7) is studied in Sec. 3. After a considerable amount of algebra, we show how this 
translates into a (strong) condition on the form of ttk- either ttk is an arbitrary ultralocal 
function, 4 or ttk must be of the form 



(d-1) 



7T K = ±2^ ^-j-J-yJW) + Cl (t)R + C 2 (t) , (8) 

where c\(t) and C2(t) are arbitrary functions of time. Comparing with (4) and (5) above, we see 
that this is of the same form as the GR result, except for the more general time-dependence. 

3 The Poisson bracket on the phase space (h i j,ir' l ' J ) is defined to be 

r , tv, f,i / SA 6B SA SB \ 

{A,B}= d d x — . . - - . . . 6) 

4 By ultralocal, we mean that ttk can depend on t and the phase space variables hij, tt^ , but not on their 
spatial gradients. 



3 



The second constraint on the allowed form for tth and tck comes from demanding that 
the momentum constraints Hi are preserved under time evolution: 

#i~0. (9) 

This condition is studied in Sec. 4. After some algebra, this translates into a (strong) condition 
on the form of the physical Hamiltonian density tth'- in the ultra-local case, 7r# and hk are 
related by a differential equation — see (74); in the case where ttk takes the square-root 
form (8), tth is also constrained to take a square-root form, and the functions ci, C2 are 
related to one another by C2 = 2A|ci| d ^ <i-1 ^. Specifically, in this latter case the result is 

2 

7l H = -U)(t)lTu , 7l K = -71^ , (10) 

where 

tt w = ±y/d{d - l)yV(2) ± S-W*) + 2Ae<M*) , (11) 

with u(t) = ^-j- log |ci(t)|. With the replacements co — > ft, A — > A, this becomes identical to 
the GR result given by (4) and (5). 

To summarize, the combined requirements of closure of the constraints and consistency 
under time evolution together single out GR as the unique effective field theory of the graviton 
degrees of freedom. We conclude in Sec. 5 with a discussion of the possible implications of 
these results for our understanding of space-time symmetry in gravitational theory. 

2 Spatially-Covariant, Transverse, Traceless Gravitons 

In this section, we construct the spatially covariant effective action of the transverse, traceless 
graviton polarizations. We describe these degrees of freedom using a unit-determinant metric 
hij subject to momentum constraints Hi. 

2.1 The dynamics of GR 

In terms of the ADM [13] variables, 

ds 2 = gflu dx^dx u = -N 2 dt 2 + (dx* + N { dt) (dx j + N j dt) , (12) 
the action of GR in canonical form is 

S= f dtd d x (hij^ - N"H^\ , (13) 
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where d is the number of spatial dimensions, 7r iJ is the canonical momentum conjugate to h^, 
and A^ = N. The lapse scalar A" appears in (13) as a Lagrange multiplier which enforces the 
vanishing of the Hamiltonian constraint 

U = 4= (hach d - j^—,h ab h cd \ n ab n cd - Vh (R - 2A) , (14) 



where h is the determinant of hij, R is the Ricci scalar of hij, and A is the cosmological con- 
stant. Similarly, the components of the shift vector N 1 appear in (13) as Lagrange multipliers 
which enforce the vanishing of the momentum constraints 

Hi = -2h l3 V k ir jk , (15) 

where the symbol V< denotes the covariant derivative with respect to h^. General relativity 
is thus properly understood as a constrained field theory. 5 Up to a sign, the Hamiltonian 
constraint (14) fixes the trace of the momentum tensor (/iy-7r y ), thereby eliminating the scalar 
polarization of the graviton. The momentum constraints fix the divergence of the momentum 
tensor (V/7r tJ '), eliminating the longitudinal polarizations of the graviton. Accordingly, the 
graviton of GR has only transverse, traceless polarizations. 

As discussed in [4], the trace of the momentum tensor is conjugate to the conformal part 
of the spatial metric, so the Hamiltonian constraint renders the conformal factor essentially 
non-dynamical. To see this, decompose the spatial metric into a conformal part e n and a 
unit-determinant part hij, 

hij = e% i3 . (16) 

Inserting this decomposition into the first term of the action (13) yields 

hij7T ij = (l (//,,-") + hij (e V J ) . (17) 

Inspecting the coefficients of Cl and hy in this equation allows us to extend the conformal 
decomposition to the momentum tensor. First, it is obvious that the conformal mode Q is 
conjugate to the trace of the momentum tensor. Second, since hij is traceless, the momen- 
tum conjugate to the unit-determinant metric is proportional to the traceless part of the 
momentum tensor. 6 We therefore define 

n n ee h^i ■ 7f« ee e n - -M 3 K b n ab 



''" ' " V" 'I 

5 See [14, 15] for pedagogical treatments of the general theory of constrained systems. 
iij follows from the identity Sh = h ■ h^Shij 



6 The tracelessness of hij follows from the identity Sh = h ■ h^Shu and the condition h = 1. 
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where 7r Jjf is manifestly traceless. The momentum tensor it^ admits the trace decomposition 

7r« = e~ n W ij + -h ij 7r n . (19) 
a 

Meanwhile, (17) simplifies to 

h i3 7i ij = (inn + , (20) 

so the spatial metric phase space (ft.^-, 7T V ) splits naturally into a conformal part (Q, ttq) and 
a unit-determinant part (fry, 7r y ). Schematically, 

(/iy,7rtf)->(fi,7rn), (kj,^). (21) 

Having performed this phase space decomposition, we will now examine how the Hamiltonian 
constraint Ho fixes the magnitude of irn and the momentum constraints Hi fix the divergence 
of 7f JJ . Substituting the trace decomposition (19) into the Hamiltonian constraint (14) yields 



(22) 



Similarly, the momentum constraints (15) can be written 



Hi = -2~h i3 V k ^ k - ^Viim ■ (23) 

Up to a sign, the Hamiltonian constraint Ho forces tiq to be a function of the remaining 
phase space variables. Since the momentum conjugate to Q is almost completely fixed by 
the Hamiltonian constraint, Q does not correspond to a true, propagating physical degree of 
freedom. At the same time, the H^s constrain the longitudinal part of 7T 1 - 7 ' to be proportional 
to the gradient of 7Tq. The only momentum variables which can be freely specified as initial 
data are the transverse components of the traceless momentum tensor tt^ . 



2.2 Our set-up 

To relax Lorentz covariance to spatial covariance, we drop the Hamiltonian constraint Ho- 
To avoid a scalar graviton polarization, we also drop the conformal part of the phase space, 
(Q,7in). We therefore choose our phase space to consist of a unit-determinant metric and 
a traceless conjugate momentum 7T V ' subject to momentum constraints Hi that enforce spatial 
covariance. The most general canonical action of this kind takes the form 

S = [ dtd d x — ix h — N^) , (24) 
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where n H is the physical Hamiltonian density. Note that has (d+2)(d— l)/2 components, 
so the spatial dimension must obey d > 2. 

In the action (24), the Lagrange multipliers iV* enforce the vanishing of the momentum 
constraints Hi- To determine the form of the Hi's, we will adapt the T^'s from (23) to the 
reduced phase space (hij,^). This ensures that our action is general enough to represent 
all spatially covariant theories of a spatial metric hij for which the conformal mode happens 
to be non-dynamical. In terms of the covariant derivative Vj with respect to the metric hij, 
equation (23) becomes 

Hi = -2hijV k f jk - (d + 2) ■ h l3 n jk V k n - 2 - ViTTn . (25) 

To simplify the form of the corresponding Hi's, we make use of the fact that we are treating 
gravity as an effective field theory to eliminate Q as an independent function. 

To treat gravity as the effective field theory of hij and 7r iJ , we must choose a background 
about which to define the action. We wish our background solution to respect the global spatial 
symmetries of the theory, so we look for maximally symmetric background field configurations 
h®j and tTq. As is well known from FRW cosmology, maximally symmetric spatial metrics can 
be flat, open, or closed. However, the fixed conformal factor of a unit determinant metric rules 
out the open and closed possibilities, so the unique maximally symmetric unit-determinant 
spatial metric is 

h% = . (26) 

Maximally symmetric rank-2 tensors living in a maximally symmetric space must either vanish 
identically or be proportional to the metric itself. The momentum tensor is traceless, so it 
cannot be proportional to the metric. It must therefore vanish on the maximally symmetric 
background, 

*? = • (27) 

Now consider the behavior of the two scalar functions Q and 7Tq on this background. To 
be consistent with spatial homogeneity and isotropy, the background scalars f2 and 7r^ can 
depend only on time, 

Q (t), vr°(t). (28) 

Now consider deformations from the maximally symmetric background. Up to this point, 
we do not have an independent definition of time, so we can use our coordinate freedom to 



7 



choose one of the scalar functions, or 7r^, to serve as a time coordinate. Choosing ttq as the 
definition of time corresponds to "constant mean curvature" gauge, while choosing Q as the 
definition of time corresponds to using redshift to define "look-back time." Either approach is 
consistent so long as the chosen function evolves monotonically. On generic backgrounds, such 
a definition of time can be applied locally, but on cosmological backgrounds such a definition 
can be applied globally. 

In this paper, we choose Q to serve as our time coordinate; inverting this relation allows 
us to write Q as a function of time. Since Q is a function of time, spatial gradients of Q vanish. 
Applying this observation to equation (25), we define 

Hi = -2/^V fc ^ - Vi7i K , (29) 

where 7Tr- is an arbitrary function of time t, the phase space variable and fr 13 ', and the 
spatial gradient operator <9j. This completes our specification of the action (24), which we 
reproduce here for convenience: 

S = J dtd d x (jiiftv -Ti H - JSMiij • (30) 

The action (30) depends on two unspecified scalar functions: the physical Hamiltonian den- 
sity nu and the momentum constraint density ttk- In the enlarged phase space (hij,^ 3 ), 
the function 7i"x would be interpreted as proportional to tiq. As we mentioned before, the 
advantage of our approach is that we can remain agnostic about how the conformal mode 
is constrained. By construction, gravitons in this theory lack a scalar polarization, and are 
therefore automatically traceless. Assuming that the momentum constraint algebra closes 
properly (see Sec. 3) and that the equations of motion preserve the constraints (see Sec. 4), 
the gravitons of this theory are guaranteed to have the correct number of polarizations. 

2.3 Building tth and tt^ 

The action (30) depends on two scalar functions tt^ and 7i"x which we will now build in 
full generality. The most general scalars consistent with spatial covariance can depend ex- 
plicitly on time t, and can be built out of three objects with spatial indices, namely the 
unit-determinant metric h^, the traceless momentum tt^ , and the spatial gradient operator 
<9j. For convenience, the metric, momentum, and gradient can be organized into the four 
spatially-covariant tensorial objects 

hij , 7T V , Rabcd , Vj , (31) 
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where R a bcd is the Riemann tensor of h^. So far, ir H and ir K are unspecified scalar functions 
of time t and the four tensorial building blocks. To enumerate all the distinct scalars on which 
tch and ttk will be allowed to depend, we again employ the logic of effective field theory. In 
determining the relevance of a given scalar, each momentum tensor is accompanied by 
a temporal expansion parameter e t , while each covariant derivative V» is accompanied by a 
spatial expansion parameter e s . For the effective field theory approach to be well-defined, we 
assume that < Et < 1 and < e s < 1. The product of expansion parameters associated 
with a given scalar will be termed its relevance coefficient. 

At zeroth order in e s , the possible scalars are ultralocal functions of the phase space 
variables, i.e., functions constructed solely out of contractions of and 7f y without any 
spatial gradients [8]. The distinct fully connected ultralocal scalars can be thought of as 
chains of momentum tensors contracted into rings. Define the scalar <p(n) to be the fully 
connected contraction of n factors of 7f* J , i.e., 

(j){n) = hij^ ij (n) , U ij (n + 1) = r a ~h ab U b3 {n) , IF(0) = h ij . (32) 

Since 0(0) = d and 0(1) = 0, any ultralocal scalar will be some function of the 0(n)'s with 
n > 2. In terms of our expansion parameters, 

<f>(n) ~e?-e° s , (33) 

so the relevance coefficient of <fr(n) is simply e™. 

As it turns out, there are no scalars of odd order in e s . This is because the spatial 
gradient <9j is the only object in the theory with an odd number of indices, so any quantity of 
odd order in gradients will always have at least one free index. 

Each scalar of second order in e s is built by contracting factors of hij and 7r 1 ^ against 
one of the following three tensors, 

Rabcd, v a v b r\ v r r b v s ri. (34) 

The leading two-derivative scalar is the Ricci scalar, given by R = h ac h bd R a b c d- in terms of 
our expansion parameters, 

R ~ e 2 s , (35) 

so the relevance coefficient of R is e 2 s . All remaining scalars of second order in e s are suppressed 
by additional factors of e t . The Ricci scalar thus represents the leading local correction to 
ultralocal physics. 7 

7 Any other scalar at the same order in e s is necessarily higher-order in e t . For instance, TT ac TT bd R a bcd ~ ^t e l- 
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To implement a consistent effective truncation of the general action, we introduce a 
truncation parameter 5 > and exclude all scalars with a relevance coefficient less than or 
equal to 5. In this paper, we wish to examine the leading local correction to ultralocal physics, 
so we assume e s < e t . To ensure that the leading local scalar (i.e., R) is not excluded, we 
choose 5 < £ 2 S . To exclude local scalars of higher order in time derivatives, we choose 5 > s 2 s e t - 
To exclude local scalars of higher order in spatial derivatives, we choose 5 > 

To determine which ultralocal scalars are not excluded, let iV momentarily denote the 
largest integer for which ef > 5. The truncated action is then allowed to depend on the 
ultralocal scalars 0(2), 0(3), . . . , 0(iV). In principle, however, N could be arbitrarily large, so 
the action could depend on any given ultralocal scalar. In practice, therefore, we will not 
restrict the possible dependence of the action on ultralocal scalars. 

Following truncation, tch and hk could depend on any of the ultralocal scalars, but can 
depend on spatial gradients only through R. In other words, 

n H = ir H (t,(fr(n),R) , it K = n K (t, (fr(n), R) , (36) 

where 0(n) stands collectively for 0(2), 0(3), . . . 

2.4 Example: Spatially Covariant General Relativity 

As an example of our formalism, let us review how to cast the space-time covariant action (13) 
of GR into the spatially covariant form (30) [4]. Following (1), we split the metric into a 
conformal part Q and a unit-determinant part hij. About cosmological backgrounds, we take 
the conformal factor Q to be our measure of time. Formally, we can use a Lagrange multiplier 
■0 to enforce a gauge-fixing constraint x = — Q(t) = 0, so that the new gauge-fixed action is 

S' = J &td d x (ir ij hij - N^n^ - ipx) ■ (37) 

The constraints "H an d X are second class, 8 so they can be solved explicitly to yield expressions 
for Q and its conjugate momentum in terms of t, h^, n 1 ^ and spatial derivatives. Substituting 
these expressions back into the gauge-fixed action yields an action of the form (30), with 

2 

■k h = -0(f)7rn , ir K = -n n , (38) 

and 

vr n = ±y/d{d - 1)^0(2) - R~^rm) + 2A e ^W . (39) 

8 Second class constraints are constraints which do not commute under the action of the Poisson bracket. 
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We verified in [4] that GR in spatially-covariant gauge is consistent as a constrained field 
theory. In the next two sections, we will show that it is essentially the unique effective field 
theory emerging from the action (30). 

3 Closure of the Constraint Algebra 

To generate a consistent gauge symmetry, the momentum constraints Hi given by (29) must 
be first-class under the action of the Poisson bracket, 

{Hi(x),H a (y)}~0. (40) 

The symbol ~ denotes weak equality, which means equality after the imposition of all con- 
straints. The condition (40) means that the algebra of the constraints is closed, in the sense 
that the Poisson bracket between any two constraints yields an expression proportional to the 
constraints. 

Satisfying the weak equality (40) will impose restrictions on the form of ttk- As we will 
see, there are only two solutions to this condition: either ttk must be an ultralocal function of 
t and the phase space variables, or ttk must take the GR form (38). However, the constraint 
redundancy implicit in the notion of weak equality makes it impossible to solve (40) directly for 
7Tk- To proceed, we must first derive a strong equation, one that holds identically irrespective 
of any constraints. 

3.1 Deriving A Strong Equation 

We defer the computation of the Poisson bracket {Hi(x), 'H a (y)} to Appendix A. The result 
is (see (152) and (153)) 

{fii{x), H a {y)} = H a (x)d x iS 3 (x -y)- U t .(y)d y «5 3 (x - y) 

+ d xi (-l k (x)d xk d x «5 3 (x - y)) - d ya (-I k (y)d yk d yi 5 3 (x - y)) , (41) 

where 
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In the expression for Z&, the symbol 5^ b denotes the traceless, symmetric Kronecker matrix, 

S = \$A + - \^h ab . (43) 
Since "Hj ~ by assumption, (41) agrees with condition (40) if and only if 

Tk ~ . (44) 

The weak equality in (44) is the necessary and sufficient condition for the momentum con- 
straints Hi to generate a consistent first class constraint algebra [4]. To solve (44), we will 
first eliminate the constraint redundancy, and thereby promote the weak equality to a strong 
equation. 

Recall that the momentum constraints % (30) fix the longitudinal part of the momentum 
tensor, Vj7f y ', in terms of the gradient of a scalar function, • In particular, the momentum 
constraints enforce the weak equality 

Vj-7f iJ ' ~ --VV* . (45) 

We therefore search the definition of Z& in equation (42) for every instance of Vj7f ij ", i.e., 
for every instance in which a derivative operator 1) acts on a momentum tensor and 2) is 
contracted against one of the indices of that tensor. Next, we perform the substitution 

->• --VVjf . (46) 

After removing every instance of Vj7p J ' in this manner, we can write Zfc ~ Z&, where 

Z fc = — ^- > m tt j V i n(m - 2) ak + — ^- > mll(m - l)^™^-— 




m II(m - 2) ka \vtck- -i^r^k > , ^^r? — 7 0( m 




^m=3 




d(f)(m) 

It follows immediately that Z& ~ 0, but we can do better. Since X k does not contain any 
factors of Vj7r^ or Hi, the weak equality Zj. ~ implies the strong equation 

Z fc = 0. (48) 

This is the desired result. This strong equation is a necessary and sufficient condition for the 
constraints to generate a consistent first-class algebra, and it contains no constraint redun- 
dancy. In what follows, we will solve this equation to determine the allowed form of ttk- 
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3.2 Solving The Strong Equation 

The simplest way to solve the strong equation (48) is for ttk to be an ultralocal function. 
Ultralocal theories of gravity have interesting applications near cosmological singularities [16, 
17] and for studying cosmological perturbations on super-horizon scales [18]. By ultralocal, we 
mean that ttk depends on t and the phase space variables hij and n 1 ^ , but not on their spatial 
gradients. In particular, dnx/dR = in this case, and X k vanishes identically. Therefore, any 
ultralocal itk = 7r^(t,0(n)) satisfies (48). 

However, in the ultralocal limit, each point in space evolves like a separate universe, so 
considerations of causality demand that we allow for local coupling between adjacent points 
in space. 9 To understand the propagation of local graviton degrees of freedom, it is of course 
necessary to go beyond the ultralocal limit. Consistent with our earlier assumptions, we allow 
txk to depend on spatial gradients only through the Ricci scalar. 

In the local case, determining the form of ttk is equivalent to solving the strong equa- 
tion (48) assuming dirx/dR ^ 0. The solution can be obtained in a straightforward manner 
analogous to the technique known as separation of variables: if one side of an equation de- 
pends on a variable X and the other does not, then in fact neither side depends on X. The 
left-hand-side of (48) is the vector quantity Xf., which depends non-trivially on an infinite 
number of distinct phase space scalars and vectors; the right-hand-side of (48) is 0, which 
does not depend on any phase space quantities. Keep in mind that equation (48) cannot be 
solved by positing additional restrictions on the phase space scalars and vectors, because this 
would reduce the number of physical degrees of freedom. Eliminating the dependence of Xj. 
on a given phase space quantity therefore imposes a corresponding restriction on the form of 
ttk- The most general ttk consistent with all such restrictions is the general solution to (48). 
In practice, it is easiest to solve the strong equation iteratively by 

1) picking a phase space quantity which appears in only a few terms of X/~; 

2) determining the corresponding restriction on ttk', 

3) substituting the restriction back into the strong equation (48), thereby simplifying it; 

4) repeating until the strong equation is solved. 

To quickly simplify (48), consider the dependence of X k on the phase space vectors 

n ja VjU{m - 2) ak , m > 3 . (49) 
9 For more on ultralocal gravity and the separate universe approximation, see [19]. 
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As is clear from equation (47), each such vector enters with a scalar coefficient 

OR d^m) ■ ^ 

To eliminate the dependence of Xk on the vectors tt^V jH(m — 2) a k without generating new 
constraints on phase space, each of these scalar coefficients must be made to vanish. Since 
Otik/QR 7^ 0, it must be the case that 

d-K K 



dcj)(m) 

In other words, the most general form for hk is 



0, m>3. (51) 



n K = n K (t, (f)(2), R), (52) 



in which case the strong equation (48) reduces to 



o ,' dTTK d 2 7T K _ dn K d 2 n K . ^ 

- i 9 r em 2 dm or* Wj4>[) 



OR 8(f)(2) dR 8(f)(2) OR 2 

8ti k ( dn K \ 2 2(d—i) d 2 ir K 
~8R \d<p{2)) + d 3(f)(2) dR 

( 8tx k (8it K V 2(d-l)d 2 7r K 
^\d^(2)\dR) + d ~d¥ 

For this equation to hold without introducing new constraints on the phase space variables, the 
coefficients of fr^V j(f>(2), 7r k 3 VjR, V^0(2) and V^i? must vanish independently. Specifically, 
the coefficients of the last two lines imply 



I V fc 0(2) 

V fcJ R. (53) 



d 2 TT K -d ( &7TK \ 2 8lT K 

d<j>{2) dR ~ 2(d - 1) \d(f>(2) J OR ° ' 

d 2 n K —d 8tt k ( diir^ ' 



1 K 



dE? 2(d - 1) d<f)(2) V dR 
Meanwhile, the coefficients of the first and third line combine to give 

d 2 n K —d ( 8tt k \ 3 

8(f)(2) 2 = 2(d-l) \d(f)(2) ) ' 



(55) 



(56) 



The fourth condition (say vanishing of the second line) follows automatically from these. To 
solve these three restrictions on n K , first write (56) as 

8 / 8n K \~ 2 d 



8(f)(2) \d(f)(2) J (d-1) 



(57) 
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Integrating this equation yields 

(H)) _2 ^( 0(2) + /(a) )' (H 

for some arbitrary function f(t,R). This is equivalent to 



^ § 1 (59) 



^< 2 > V d ^(2) + /(*,«) 
Integrating this equation yields 



n K = g{t,R) ± 2yJ^-^-y/<j>(2) + f(t,R), (60) 

where we have introduced a second arbitrary function, g(t,R). This form for ttk solves (56) 
in full generality. Using (60), condition (54) becomes simply 

which implies that g = g(t). However, since ttk appears in the ~Hi constraint through V^ttr-, 
the function g(t) drops out of the action (30) and can therefore be set to zero without loss of 
generality. Using (60) and g = 0, the final condition (55) becomes 

d 2 f 



dR 2 



0, (62) 



the solution to which is 



f(t,R)= Cl (t)R + c 2 (t), (63) 



where c\(t) and c^it) are arbitrary functions of time. The most general solution to all the 
restrictions on ttk is therefore 



(d-1) 



w K = ±2yJ ^-^0(2) + c x (t)R + c 2 (t) . (64) 

This n K has the same dependence on the phase space variables as the n K of GR — see (38) 
and (39). The only difference lies in the explicit time dependence. The ttk of GR depends 
on only one arbitrary function, Q(t), while the itk we just derived depends on two arbitrary 
functions, c\(t) and C2(t). These will be further constrained in the next Section by demanding 
that the constraints are consistently satisfied under time evolution, leaving us with GR as the 
only possibility. 
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4 Evolution of the Constraints 



For the constraints Hi to be consistent with the equations of motion, they must be preserved 
under time evolution, i.e., 

-Hi-0. (65) 

By promoting this weak equality to a strong equation, we will derive a restriction on the 
allowed form of the physical Hamiltonian density 7r#. 
The Hi generators evolve according to 

~ dV(x) 

H i (x) = - t ^- + {H i (x),H}, (66) 

where H = J d d x (tth + N l l-ii) is the Hamiltonian. Since the constraints are by assumption 
first class, {Hi(x), 7ii(y)} ~ 0, and therefore 

~ dW(x) 

Hi(x)~-?^ + {Hi(x),IL}, (67) 

where 

n = d d x7c H . (68) 



Because II coincides with H when the constraints are satisfied, i.e. II ~ H, and because 
it contains no constraint redundancy, II is called the physical Hamiltonian. Noting that %i 
depends explicitly on time only through tik, (67) reduces to 

iUt?) ~ ~ d i {^T-) + Wiix), n} . (69) 

To proceed, we must calculate {U,'Hi(x)}. The derivation is rather technical, and we leave 
the details to Appendix B. Substituting the resulting expression (165) in (69), we obtain 

<ti ( \ a ( , d7r K (x) 2dir H ~ ^ d-K H , f v 2(d - 1) ~ ~ b d7i H 

— fidx) ~ di its H tt 1 -^- — ; — —R - > n - , . . din) — - — -V fe V — — 

K J \ dt d 8R ^ d(t){ny y 1 d OR 

r\ OO j-» r\ OO r\ 

H —R > "^7777 — rllfm— lj6c<3^ —R > n ni , , Hn - ljjfcfT 

OR ^ d<P(m) ] dR ^ d<t>{n) ,] 1 



be 



<90(m) lbc > k \dR) 



m=2 

r> / OO . 
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On the right hand side, everything appears beneath a single gradient operator. The consis- 
tency condition Irii ~ therefore implies that the quantity on which the gradient acts must be 
weakly equivalent to some function of time; call it f(t). However, we are free to redefine the 
quantity ith — f(t) to be ti h ; this redefinition changes the action only by a total derivative, and 
leaves the equation of motion invariant. This eliminates the arbitrary function f(t), allowing 
us to write 

dir K (x) 2dn H ~ ^ dn H 2(d - 1) ~ ~ b dn H 

~ H o ; — — R — > n - , . . (bin) ; V;,V — — 

dt d dR ^ d<f){ny K ! d dR 

+ t^E^n (m - 1)J% - E^n (m - ^ (H) 

To promote this weak equality to a strong equation, the key observation is that gradients act 
on loose momentum indices only in very the last term, thus all of the constraint ambiguity 
arises from this term. Expanding the last term, applying the identity ViVjM ab = VjViM ab + 
R a cij M cb + R b cij M ac , and performing the substitution — > — |VV^, we obtain the strong 
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equation 

dir K (x) 2dn H ~ ^ dn H , f . 2(d - I) ~ ~ b dn H 

= tth H ^ ; — —R — > n———(p(n) V;,V — — 

dt d OR ^ d(f){ny K ' d dR 

d dR ^ d<j)(ny K ' d 8R ^ <90(m) rv ' 

- sir (fi "lR n( " - 2 »-) vt - 2 lr (*~) ^ (m) 

+ i H H (VV) (f »^n (n - 3).) - if! H ^ (£ ,^n (n - 3).) 

(72) 

We next import our knowledge of hk from Sec. 3 to simplify (72). Recall that hk can either 
take the ultra-local form tt^ = t^k^i 0( ri ))j or the square-root form (64). We will treat each 
case separately. 

4.1 Ultra- local case 

Since ttk = vr^(t, 0(n)), and in particular di^x/dR = 0, in this case the strong equation 
simplifies to 

dn K 2dn H ~ ^ d-K H \ , dir H ~ jk ^ dn K 

= ir H + — — R— > n———d>[n) H ^i?- 7 > m— - — -Il(m - l) ik 

dt ddR 4^ ; dR ^ 9 d0(m) 1 

ld7c H ~s^ dit K 2{d- 1) ~ ~ k dn H 

; — —R / m— — — -0(m — 1) : VfcV — — 

d 8R ^ d<Hmy K ' d dR 

m=3 

+ 3 r ™#^(™ - 1) V*V* (^)-T -#^n(m - l^V, V, ( . (73) 
rfZ-, d<P(mf K } \dR) ^ <90(m) 1 ; 3 k \dRj K ' 
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As it turns out, equation (73) can only be solved if n H is ultralocal. To prove this, we use 
the same techniques involved in the derivation of equation (64). We defer the details of the 
analysis to Appendix C. 

Since n H must be ultralocal, d^n/dR = and the strong equation (73) becomes 

° = ^ + ^f-E-^H- ( 74 ) 

This is the analogue of the conformal renormalization group (RG) equation discovered in [4]. 
It implies that the time-independent part of ttk can be arbitrary, whereas the time-dependent 
part is fixed by 7r#. To summarize, the ultra-local case allows for two arbitrary ultralocal 
functions of the phase-space variables: a possibly time-dependent tth = 7i\ff(£, <fi( n )) an d a 
time- independent component of ttk- This completes our analysis of the ultralocal case. 

4.2 Local case 

We turn now to the local case, in which 7Tk takes the square-root form given by (64): 

n K = ±2^j Mzi) ^0(2) + Cl (t)R + c 2 (t) . (75) 

We of course assume that c x ^ 0, so that dnx/dR ^ 0, for otherwise this would reduce to the 
ultra- local form. Since tck is a function only of t, R, and 0(2), the strong equation (72) can 
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be written as 

dir K (x) ^ dn H 2dn H ~ ldir K ~^ dn H dn K dir H „ ~ t 

= tth H - > n———(b(n) — — i2 + - — > n - , . . 0(ra - 1) - 2 — — - - . . ir bc R 

dt ^ d(f){ny K ' d dR d dR ^ d<j){ny K ' OR 50(2) 

- ^ (g - ^ (H) - O 

+ iff M N (|^n ( „ - a).) - ^ H *. (H) 

- ((^-) g-fe*- - + ^ (H)) 

+ If (^) (g "l^ n( " " 2) ^) + lr" M (g "l^> n( " " 

The phase space scalars 



(76) 



« db V c V b n(n-3) ad , n>4 (77) 

appear only in the final summand of (76). To remove these scalars from the right-hand-side 
of (76) without introducing new constraints on phase space, their respective coefficients 

dir K dir H 

— ^n— — — , n > 4 (78) 
dR d(j){n) 

must vanish. Since diix/dR ^ 0, we must have 

-Jtt^t = °) n>4, (79) 
o<p[n) 
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so n H can depend only on t, R, 0(2), and 0(3). The strong equation (76) then becomes 

dn K (x) ^ 8-kh v 2d-K H ~ ld-K K ~d-K H ,, n , Q dn K dir H _ ~ 6c 

= ^ + -at— 5 n ^ 0(n) ~ Sa3F* + W^sf (2) - 2 a*pf^ 

dRdct>{2) \ cd ab ) dRd<t>(3) d dR d dR W(3) 

+ 2^"Vj% ( g/L) - 2#^V,V fc Ynp^Tl(n - 2)A V a V b n K 

OR 3 \d<j>(2)J 90(2) 3 \dR) 2 dR d<j>(n) 1 ' ) 

~ ^ (--) ^ (H) + 3 iH (v-) H " 4 1 (*~) ^ (H) 

-4t^-(M (H)) +3 f ^HM) -H^^H)) 

To remove the scalar ( V b 7f ac ) ( V c 7r a b) from the right-hand-side of this equation without intro- 
ducing new constraints on phase space, its coefficient 

d% K dir H 
dR d<t>{Z) 

must vanish. Since dn^/dR ^ 0, we must have 

dir H 



1 



90(3) 



0, 82 



so txh can depend only on t, R, and 0(2). Using this fact, and making use of the explicit form 
of 7Tk in all partial derivatives of ttk, the strong equation (80) becomes 

- - M + 2 (If -m) *>* 
-"**{%-«mY (83) 



By defining 



s ig-«m (84) 
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- Ci 



and using the chain rule, we can write the strong equation (83) as 

d ~ dir H d dn H ~ 1 ~ ~ ditu 

= W^if" 1 "' + Cl +C2 ~ m) m ¥^) WK ~ ~dP R ¥^Tf' - " kV ' v wf 

- ^ t (Vi g^) VrfW - * Jg**W - 2 g ^f^V^M . (85) 
To eliminate the dependence of the right-hand-side on the scalars 

Tt^VjVkR, n jk VjV k (j)(n), n>2 (86) 

without imposing any new constraints on phase space, the respective coefficients of the scalars 
must vanish, 

^£ = 0, tt^t=0, n>2. (87) 
OR d<j>{n) ~ y ' 

It follows that / = f(t), which from equation (84) implies that 

dn H d7T H 

-3R =Cl (m) +m - (88) 

Substituting this result into the strong equation (85) and simplifying a total derivative yields 



d ~ dun /riN d ( diiu . \ ~ 1 



39) 



By defining 



g (t,R,4>(2))=7r K -^- (90) 



and using the chain rule, we can write the strong equation (89) as 



2(^1)™ + ClR + C2 ~ dm m ¥^~f K ~[ Cl dm + m ) Jd^T) 



V <w 50(n)y or t^2 d ^ 

(91) 

To eliminate the dependence of the right-hand-side on the scalars 

V a V a £, V a V a 0(n), n>2 (92) 
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without imposing any new constraints on phase space, the respective coefficients of the scalars 
must vanish, 

= , = , n > 2 . (93) 

OR d<j>{n) ~ y ' 

It follows that g = g(t), which from equation (90) implies that 
or 



±9(t)-.l _, - , (95) 



^< 2 > 2 V^-i)^ (2)+Cl(i) ^ +C2(i) 

Integrating once yields 



7c H = h(t, R)±g{t) x j y / 0(2) + c 1 (t) J R + c 2 (t) , (96) 

where R) is an arbitrary function. Taking a derivative with respect to R, we find that 



Sfl Sfl l9W 2\/(rf-l)./ 0(2)+Ci(t)i j + C2(t) ' 



Comparing with equation (88) yields 

dh 



OR f(t) ' <* 



so 



h(t,R) = f(t)R + k{t) , (99) 

where k(t) is an arbitrary function. The physical Hamiltonian density tth can therefore be 
written as 



k H = k(t) + f{t)R±g(t) x j y / 0(2) + c 1 (t) J R + c 2 (t) (100) 

However, the function k(t) contributes only a boundary term to the action and does not affect 
the equations of motion. We can therefore choose k = without loss of generality, in which 
case 



n H = f(t)R ± g(t) x j ^-yy\M 2 ) + ^R + c 2 (t) . (101) 
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With this result, the strong equation (91) becomes 

= 2?d~-l) * Kf ®** + 9{t)ci{ t )k + 9{t)^{t) t^I) + ^ + ^ + 2 *jkR jk f(t) ■ (102) 

If d > 3, eliminating the phase space scalar TTj^RP from the right-hand-side of the strong 
equation without generating new constraints requires / = 0. If d = 2, then R = yhR in 
equation (101) is a total derivative, and once again we can choose / = 0. Since / = either 
way, Tin is 



n H = ±0(^/^-^0(2) + Cl (t)R + c 2 (t), (103) 

and the strong equation becomes 

= c 2 + c 2 —^--g(t)+(c l + c 1 g(t)}R. (104) 

To eliminate the dependence of the right-hand-side of the strong equation on R without 
reducing the number of phase space degrees of freedom, its coefficient must vanish, 

c 1 +g(t)c 1 = 0. (105) 

Since c\ ^ by assumption, we conclude that: 

g{t) = --- (106) 

Cl 

Meanwhile, the remainder of (104) has the general solution 

c 2 = 2A|ci|^, (107) 

where A is a constant. 

These results are most neatly expressed in terms of a new time-dependent function: 

Wft^^-y logldM. (108) 

Writing C\ and c 2 in terms of uj allows us to define the scalar function 

tt w = ±y/d(d - 1)^0(2) ± Be^-iMt) + 2\ e Mt) , (109) 
in terms of which 7T# and hk can be written as 

2 

7T^ = -w(t)7T w , 1T K = -7T W . (110) 

d 
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Apart from the sign of the coefficient of R in ttq, these last two equations are identical to (38) 
and (39) for GR. Since the R term determines the 'gradient energy' of gravity, the sign 
can be fixed by requiring the absence of gradient instabilities. For instance, on looking at 
tensor perturbations around spatially flat FRW backgrounds, a gradient instability will only 
be absent for the choice of sign corresponding to GR. The function u(t) and the parameter A 
are identified as u(t) and A of GR, respectively. 

4.3 Summary 

Demanding the consistency of the momentum constraints Hi with the equation of motion 
yields only two possibilities for tth and ttk- First, hh and tt^ can be ultralocal functions of 
time t and the phase space variables, subject to the restriction 

In this case the physical Hamiltonian does not include spatial derivatives of the fields, so there 
is no notion of gradient energy, and no physical coupling between adjacent spatial points. Such 
theories cannot describe the local propagation of graviton degrees of freedom. 
The second possibility is to have 

7T H = -u(t)ir u , n K = ^7T W , (112) 

where u(t) is an arbitrary function of time, and 

tt w = ±y/d{d - l)yj(f>{2) ± £e(<i-iM*) + 2Ae^W . (113) 

To leading order in local scalars, consistency alone has forced the functions tth and tcx to take 
essentially the same form as in GR. The local functions tih and tt^ obey an extended version 
of the restriction (111), namely 

In terms of 7r w , this becomes the RG flow equation 

0= I 5-^ + > 'n4>(n)-^ + -R^ 1 e"^ 2 ^ . (11.-)) 




As discussed in [4], this equation encodes the invariance of the physical Hamiltonian density 
under flow through the space of conformally equivalent metrics. 
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5 Conclusion 



In this paper, we began with the spatially covariant action of the transverse, traceless graviton 
degrees of freedom. We then applied the formalism developed in [4] to determine under what 
circumstances the momentum constraints of the theory would 1) generate a first class algebra 
and 2) be preserved under time evolution. To leading order in local scalars, we found that 
consistency alone singles out general relativity as the unique effective field theory of the 
graviton degrees of freedom. 

To our knowledge, this represents an enormous advance over all previous derivations of 
general relativity from the graviton degrees of freedom, which assume Lorentz covariance at 
the outset [1-3]. Our approach relies on the weaker assumption of spatial covariance, and yet 
achieves an equally powerful result. 

In light of our result, it is plausible to interpret Lorentz symmetry in the gravitational 
sector as an accidental or emergent symmetry. Accidental symmetries arise in an effective 
field theory when all the allowable operators which violate the symmetry are confined above 
the cutoff of the theory; in this respect, it is the opposite of spontaneous symmetry breaking. 
This is exactly the case in our derivation of general relativity, where consistency forces the 
operators 0(2) and R to appear in a Lorentz covariant combination, and our cutoff excludes all 
higher order operators which might spoil the symmetry. Also of note is the fact that both the 
conformal scale factor Q(t) and the cosmo logical constant A arise in the theory as constants 
of integration with respect to phase space scalars. 

A clear direction for future research would be to generalize the derivation by allowing 
the functions tth and ttk to depend on spatial derivatives through scalars other than R. There 
are only two logical possibilities: 1) Lorentz covariance survives at all orders in perturbation 
theory, or 2) Lorentz-violating scalars can modify the behavior of the graviton. Either way, 
the results should be interesting. 

So far, our discussion has focused solely on the graviton degrees of freedom, but any 
Lorentz covariant field theory could be subjected to a similar analysis. Currently, we are 
investigating Hamiltonians which yield a consistent dynamical evolution for spatially covari- 
ant Yang-Mills fields [20]. To systematically relax the assumption of spacetime symmetry 
throughout particle physics, it will eventually be necessary to expand our gravitational frame- 
work to include coupling between the graviton and all the relevant matter fields. 
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A Computation of {Hi(x), 7i a {y)} 

The Poisson bracket of the theory is defined to be 

{A,B }s f#x( ~ SA t . s »-J*- tB V (116) 

To evaluate the variational derivatives inside the Poisson bracket, one must use the relations 

5h rs (z) dn rs {z) 

and 

**M = -h"r-5\x- Z ), ||^=0, (118) 
Sh rs {z) d 57T rs {z) 

from which follow the operator identities 

The symbol S^ b denotes the traceless, symmetric Kronecker matrix, 

€ = ^-^^ b , (120) 
while S^ b denotes the symmetric Kronecker matrix, 

C = (121) 

To compute the Poisson bracket {Hi(x), 7i a (y)}, we first use the functionals 

F = J d d x fiU , G = J d d y g a H a , (122) 

to compute 

{F,G} = J d d xd d y r(x)g a (y){'H l (x),n a (y)} . (123) 
Write the momentum constraints in the action (30) as 

Hi = J i + )C i , J i = -2h ij Vk* ik , Ki^-ViTtK- (124) 
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This allows us to write F = Fj + F K and G = Gj + G K , where 

Fj = J d d x fJi , F K = J d d x fid , 



and 



Gj = I d d yg a J a , G K = / d d y g a K a 



It follows that 

{F, G} = {Fj, Gj} + {Fj, G K } + {F K , Gj} + {F K , G K } 

Recall that Ji = —2hijV k^ k ■ It is straightforward to verify that 

ip _ / _ \ _ 9 



0h rs V 7 « 



and 



SG 



rs 



J _ r6c 



~5 r a t (2n bc V c g a ) - V a (^tt") - ^ rs V a g a , 
S b r c s {2h ab V c g a ) ■ 



5n 

Recall that /Q = — V^k- If tTk = ^xit, <f>(n), R), then 
SF K 



bh rs u 2 



P = (W)f:»^«fn(»-i)M, 



n=2 



and 



- ( d a g a ) d -^~5ZR bc + cv 6 v c (W)f|) 



28 



Combining the Fj and Gj variations into the bracket {Fj,Gj} yields 

{Fj, Gj} = 2 J d d z { (v c f ) (V^ a ) h ab K bc - (V k g a ) (v a f ) /^' fc 

+ (V fc f ) V a ^ a %7f ifc ) - (V c £ a ) V, (/V^f*) } . (132) 

After integrating by parts, using the definition Ji = — 2h i j'V , and using the identity 
(ViVj - VjV;) V a = R a bij V b , this reduces to 

{Fj, Gj} = J d d z {fj a V ig a - g a JiV a f + 2fg a n jk [k jika + R jaik ) } . (133) 

From the symmetries of the Riemann tensor 10 and the traceless momentum tensor 11 , it follows 
that ^ (R 3ika + R jaik ) = 0, so the last term in the integrand vanishes. The connection terms 
inside the remaining covariant derivatives cancel to yield 

{Fj, Gj} = j d d zfj a d i g a - J d d zg a J i d a f. (134) 
To extract the bracket {Ji(x) , J a {y)} from this result, first relabel dummy indices 

{Fj,Gj} = J d 3 xfj a d l9 a - J tfyfJidaf*. (135) 
Under the spatial derivatives in this equation, insert the identities 

g«{x) = j d 3 y5 3 (x-y)g a (y), f\y) = J d 3 x 5 3 (x - y)f (x) , (136) 

to obtain 

{Fj,Gj} = J d 3 xd 3 yf\x)g a (y) (j a (x)d xi 5 3 (x - y) - Ji(y)dy6 3 (x - y)) . (137) 
This yields the distributional identity 

{Ji(x), Ja(y)} = Ja(x)d x ,5 3 (x -y)- J t (y)d y ^ 3 (x - y) . (138) 
Combining the Fj and Gk variations into the bracket {Fj, Gk} yields 

{Fj,G K } = - J d d z(dj%d a g a )f^n-^<t>(n) d d z (d.f )(^)^p 

- j d d zf\d a g a )^^- ) d i <j>{n) + 2 J d d z{d a g a ){R lk V k r 



OR 



-2jd d z (v fc f ) V,V fc (W)f|) + \ j d d z )V C V C (W)f|) ■ 

(139) 



Rabcd — Rcdabi Rabcd — Rbacd Ra 



bdc • 
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Using the identities VjVfc/ 1 = V^Vj/ 4 + Rikf 1 and 2V k Rik = Vj-R, integrating by parts, 
simplifying a total derivative of ttk, integrating by parts again, and expanding yields 



{F J} G K } = -J d d zf(d a g a )d t 7r K + 2{d d 1} J d d z(d a g a ) (v fc («9 4 f )) V A 

+ 2 ^/d^^(w ) ) (vw } ; 



-j^ {aM) ^R + P^ in) y (140) 

and a parallel expression for {Fk,Gj}. The sum of the two brackets reduces to 
{Fj,G K } + [F K ,Gj] = - f d d z (d„g")fd,* K + f d d z (dJ^g'd^K 



+ ^/d^)(v l( ^)v'§S 



d 

After integrating by parts, expanding, and simplifying, we obtain 
{Fj,G k } + {F k ,Gj}= [ tfzpKAcT- [ d d zg a K l d a f l 



where 



J d d z ftf) (V k (d a g a )) V fc || . (141) 



+ J d d z {d i f i )[d k d a g a )M h - J d d z (d a g a )(d k d i f)M k , (142) 



2(d - 1) ~ dir K . . 

M^-^V^. (143) 

To extract the bracket {Ji(x), JC a (y)} + {K.i(x), <?a{y)}, integrate by parts, relabel dummy 
indices, and insert the identities in equation (136) to yield 

{Fj, G K } + {F K , Gj} = J d 3 xd 3 yf(x)g a (y) (lC a (x)d x ,5 3 (x - y) - K^d^S^x - y)) 

+ J d 3 xd 3 yf{x)g a {y)d xi (-M k (x)d xk d xa 5 3 {x - y)) 

- J d 3 xd 3 yf i (x)g a (y)d r (-M k (y)d yk d yi 5 3 (x - y)) . (144) 

It is clear that 

{Ji(x), K a {y)} + {)d(x),J a (y)} = lC a (x)d x >5 3 (x - y) - K i (y)d y a8 3 (x - y) 

+ d x i(-M k (x)d xk d xa 8 3 (x-y)) 
- d y a (-M k (y)d yk d yi 5 3 (x - y)) . (145) 
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Combining the F K and Gk variations into the bracket {F K , Gk} yields 

9tt k 
dR 

9tik 
OR 



{F K ,G K } = J d d z(d a g a ) ^ m ^^n(m-l) 6c ^ VW fc 

- J d d z (dif) (^™J^n(m - l)J>tj V J V fc [(dag* 



Integrating by parts, expanding, and simplifying yields 

dR 



{F K ,G K } = J d d z(Vd t h(d a g a ) \^m^^U(m - 1)J%) 



Define 



Then 



(146) 



dR 1 9 d 0( m ) 

\m=2 v ' 



t™^*"-^)^ (148) 



{Fki G k } = J d d z (d i f)(d k d a g a )Af k - J d d z (d k d i f)(d a g a )Af k (149) 

To extract the bracket {/Q(x), JC a (y)}, integrate by parts, relabel dummy indices, and insert 
the identities in equation (136) to obtain 



- [ A 3 x<??yf'(x)f(y)d r (-M"{y)d v *d s ,S z {x-y)). (150) 



It follows that 



{ld{x), K a (y)} = d xi (~Ar k (x)d xk d x a5 3 (x - y)) 

- dya {-U k {y)d yk d,^\x - y)) . (151) 
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Since "H, = Ji + /Q, combining equation (138) with equations (145) and (151) yields 

{Hi(x),ii a (y)} = H a (x)d x i5 3 (x - y) - / H i {y)d y a5 3 (x - y) 
+ d xi (-l k (x)d xk d x «5 3 (x-y)) 

- d ya (-l k (y)d yk d yi S 3 (x - y)) , (152) 
where = M.k + A4- In terms of tik, 



^r^ {m ] ^r ^s~i Vi ^- (153) 

B Computation of {n, Hi(x)} 

Recall that the physical Hamiltonian is 

n = d d xn H . (154) 



If 71h = TTH(t, 0(n), R), the variational derivatives of IT are 

Recall from equation (124) that 

Hi = J i + )C i , J i = -2h i jVk ! K ik , K i = ~V l 7iK. (156) 

To compute {H,%i(x)}, we first compute the brackets {II, Ji(x)} and {H, JCi(x)} . Combining 
equations (155) and (129) yields the Poisson bracket 

+ y d ^2(v t9 »)v.V'(^). (157) 
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Integrating by parts, using the identity V a Vb9 a = V&V ' a g a + g a R a b, and integrating by parts 
again yields 

+ / d^2 9 "||m s + / d", (v„ 9 °) V t V' (||) . (158) 

Since 2V*£* = V„fl = d„fl and V„0(n) = 9„0(n), 

+ ^/d^(v^)v^(^). (15.) 

Simplifying a total derivative and integrating by parts yields 

{ n, Gj} = /d^ ^-/dS 9 ^(?^ + g^(n)) 

-/^^v>H). " 2 (160) 

It follows that 

{ n, j i} = ^ \n H -- dW R- E n 9$w*W ' ~^ VbV -qS ) ■ ( 161 ) 

Combining equations (155) and (131) yields the Poisson bracket 

{u,g k } = -j i«z(d„f) K>^_n(m - i)J%\ d jSw* 
+ / AM (|»^ n( "- 1) ' S: )^ 

- / (g"^n ( „ - 1)^) V»V< (W)f|) (162) 
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Repeated integration by parts yields 

// o CO o 

d d zg a d a -^§R? k V m — ^n(m - 1) 6 J} 



_/d^g m ^n (m -i)J^^ 
+ /d',^(^w(| B ^n(„-i)^ 



It follows that 



m=2 v ' x ' 

Combining equations (161) and (164) yields 

m <£} i « 2 ^5 d7t H , 2(d - 1) ~ ~ b dir H 



n=2 

oo 



oo r\ oo 

+ ft ( ^W k £ -^n(m - l)J% - V n^n(n 

1 OR ^ d<j)(m) J dR ^ d<f>(n) 

90(m) jfc v^y 



m=2 

oo 
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C Proof that ultralocal ttk requires ultralocal tth 



In this section, we prove that if ttk is ultralocal, then tth must be ultralocal. When ttk is 
ultralocal, the two functions hk and tth must satisfy equation (73), namely 

= tth H — ^ ; — ^ix — > re - , . . 6{n) H —R J > rn—— — -U(m — l) ik 

dt d dR 4^ d(t){nY y ' dR ^ d(j){m) 1 )jk 

1 Otth 5 ^ 07TK 2(d- 1) ~ ^ k d-K H 
- — — R > rn——— — -0(m - 1) — - — -V fc V — — 

m=3 

The right-hand-side of this equation depends on an infinite number of distinct phase space 
scalars, but the left-hand-side is 0, which does not depend on any phase space scalars. It 
follows that the right-hand-side cannot depend on any phase space scalars either. We cannot 
impose additional restrictions on the phase space scalars themselves without reducing the 
number of physical degrees of freedom, so eliminating dependence of the right-hand-side on 
a particular phase space scalar can only impose restrictions on ttk or tth- Instead of trying 
to derive all the restrictions at once, we will focus on eliminating the dependence of (166) on 
a few types of scalars at a time. Using the chain rule, the final summand of (166) can be 
written as 

m———U(m - lV fc V 7 -V fc — — = > m— — — — —U(m - lV ft V 7 Vfci? 
m=2 94>{m) 1 ; 3 \ dR ) ^ d(j>(m) OR 2 K J 3 

OO OO t\ rj2 

d<p(m) d( f>(n)dR 1 3 kn } 



^ d<f>(m) K } V 3 dR? J 

°° OO / r>2 \ 

y m iS^-M m - 1)* y ( v,-^u ) v fc 0H . 



(167) 



By inspection, the scalars 

n(m - l) jk VjV k R , n(m - l) ifc V i V fe 0(re) , m, re > 2 (168) 

appear in the strong equation (166) only inside this final summand, with scalar coefficients 

dir K d 2 ir H dir K d 2 ix H 

Tfl — — 777, — ( 1 Oc/ ) 

dcf>(m) 8R 2 ' d(j){m) d<j>(n)dR ' 
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To satisfy the strong equation (166), these coefficients must vanish. To prove that n H must be 
ultralocal, we will examine two mutually exclusive and exhaustive cases: first, the case when 
d7iK/d(j)(n) = for all n; second, the case when dnx/dip^n) ^ for some n. 



C.l d7Ti(/d(J)(n) = for all n 
One way the coefficients (169) can vanish is to have 

dir K 



d(j)(m) 



0, m>2, (170) 



in which case ttk = T^xit)- Since tck only appears in the action through a term Vi^x inside 
Hi, a purely time-dependent tt^ drops out of the action. Without loss of generality, we can 
set 7Tk(£) = 0. Using ttk = and the chain rule, the strong equation (166) becomes 

2dir H ~ ^ d-K H , ( , 

To eliminate dependence of the right-hand-side on the scalars VfcV fc i? and VfcV fc 0(n) without 
generating new constraints, their respective scalar coefficients must vanish, i.e., 

^ = 0, ^ = n>2. (172) 

OR 2 d<p(n)dR 

It follows immediately that 

dir H 



, p f(t), (173) 

where f(t) is an arbitrary function of time. Integrating once more yields 

n H = f(t)R + g(t,<j>(n)), (174) 

where g(t, <fi(n)) is an arbitrary ultralocal function. With this result, the strong equation (171) 
becomes 



= — f(t)R + g(t, 0(n)) - n ^h^ ■ ( 175 ) 
d f-^ d<t>{n) 
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To eliminate the dependence of the right-hand-side on R, its coefficient must vanish, 

^p/(*) = 0. (176) 

If d ^ 2, then /(£) = 0. From equation (174), it follows at once that tth is an ultralocal 
function, which is the desired result. 

If d — 2, then f(t) can be arbitrary, so equation (174) appears to imply that tth can 
depend on R. However, tth contributes to the action only through a term 



d d xn H = J d d xg(t,(f)(n)) + f(t) J d d xR. (177) 

Since h = 1, we can equally well write this as 

J d d xn H = J d d xg(t,0(n)) + f(t) J d d xVfiR. (178) 

When d = 2, the Gauss-Bonnet theorem tells us that vhR is a total derivative, so the latter 
integrand does not affect the equations of motion. Once again, we can choose f(t) = without 
loss of generality; equation (174) then implies that tch is an ultralocal function, which is the 
desired result. 

C.2 Case II: 87Tk/ 'd(j){n) ^ for some n 

If d7iK/d(j)(n) 7^ for some n, the vanishing of the coefficients (169) requires 

—J^- = 0, n>2. 179 

dR 2 d<j)(n)dR 

It follows immediately that 

5| = /(*) , (180) 
oR 

where fit) is an arbitrary function of time. Integrating once more yields 

n H = f(t)R + g(t,(f>(n)), (181) 

where g(t,<f){n)) is an arbitrary ultralocal function. 

If d = 2, the Gauss-Bonnet theorem tells us that R = V%R is a total derivative. Since 
7Th only appears in the action through J d d xiTH, without loss of generality we can choose 
f(t) = in equation (181). It follows immediately that tth is an ultralocal function. 
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If d > 2, applying (181) to the strong equation (166) yields 

- -f(t)R m-^-J>{m - 1) + f(t) m 7nT^ u ( m - l )okR jk 
d ^ d<P{m) ^ d(f>(m) 



:i82) 



To eliminate dependence of the right-hand-side on the phase space scalars 

II(m - l) jk R jk , m > 2 (183) 

without generating new constraints, the corresponding scalar coefficients 

f(t)m^- , m > 2 (184) 
<70(mJ 

must vanish. Since 7^ for some m, the function f(t) must vanish. Equation (181) then 
implies that tth is ultralocal. This concludes the proof that if 7Tk is ultralocal, 7Th must be 
also be ultralocal. 
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